set.seed(2027)
n <- 396
baseline_true <- rnorm(2 * n, mean = 185, sd = 20)
baseline_obs <- c(
baseline_true[1:n] + rnorm(n, mean = -10, sd = 20),
baseline_true[(n+1):(2*n)] + rnorm(n, mean = 0, sd = 20)
)
group <- rep(c("Semaglutide", "Placebo"), each = n)
treat <- ifelse(group == "Semaglutide", 1, 0)
true_effect <- 10
followup <- baseline_true + true_effect * treat + rnorm(2 * n, mean = 0, sd = 20)
change <- followup - baseline_obs
df <- tibble(group, treat, baseline_true, baseline_obs, followup, change)